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I discuss the saturation regime of QCD: the weak-coupling regime that describes large parton 
densities inside hadrons and the resulting high-energy behavior of scattering amplitudes. I briefly 
review past successes in the context of deep inelastic scattering at HERA and forward particle 
production at RHIC, and I present some recent progresses that significantly improved our under- 
standing of high-energy scattering in QCD. 
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The QCD description of hadrons in 
terms of quarks and gluons depends on the 
processes considered and on what part of 
the hadron wavefunction they are sensitive 
to. Consider a hadron moving at nearly 
the speed of light along the light cone di- 
rection x + , with momentum P + . Depending 
on their transverse momentum kx and longi- 
tudinal momentum xP + , the partons inside 
the hadron behave differently, reflecting the 
different regimes of the hadron wavefunction. 

When probing the (non-perturbative) 
soft part of the wavefunction, corresponding 
to partons with transverse momenta of the 
order of Aqcd^200 MeV, the hadron looks 
like bound state of strongly interacting par- 
tons. When probing the hard part of the 
wavefunction, corresponding to partons with 
kx 3> &.QCD and x < 1, the hadron looks like a 
dilute system of weakly interacting partons. 

The present work deals with the small— x 
part of the wavefunction and with the so- 
called saturation regime of QCD. When 
probing partons that feature kx S> Aqcd 
and the effective coupling constant 

a s log(l/x) is large, and the hadron looks like 
a dense system of weakly interacting partons, 
mainly gluons (called small— a; gluons). 

The larger kx is, the smallest x needs to 
be to enter the saturation regime. As pic- 



tured in Fig.f, this means that the separa- 
tion between the dense and dilute regimes is 
characterized by a momentum scale Q s (x), 
called the saturation scale, which increases 
as x decreases. The scattering of dilute par- 
tons (with kr ^S> Q s (x)) is described in the 
leading-twist approximation in which they 
scatter incoherently. By contrast, when the 
parton density is large (fcr~Q s (i)), partons 
scatter collectively. The saturation regime 
of QCD is the perturbative regime that de- 
scribes this collective behavior. 

To be sensitive to the small— a; physics, 
high energies are needed. Hence before the 
mid 90's, saturation was not relevant for any 
experimental measurements, contrary to soft 
QCD (for instance in hadron-hadron elastic 
scattering) or hard QCD (for instance in top 
quark production). But with the start of 
HERA and then RHIC, there has been a re- 
cent gain of interest for saturation physics as 
observables sensitive to the small— ir part of 
the hadron wavefunction could be measured. 

One of the most studied process is deep 
inelastic scattering (DIS): a photon with vir- 
tually Q 2 probes the hadron wavefunction 
and observables are sensitive to values of x as 
small as XBj—Q 2 /s. At high energy s^>Q 2 , 
both inclusive, diffractive and exclusive pro- 
cesses are probing the saturation regime. 
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Fig. 1. Diagram in the (kx,x) plane picturing the 
hadron in the different weakly-coupled regimes. The 
saturation line separates the dilute (leading-twist) 
regime from the dense (saturation) regime. 

In DIS, the hadronic scattering can be 
viewed as that of a colorless qq pair (or 
dipole) of transverse size r ~ 1/Q off the 
hadron. The evolution of the dipole scat- 
tering amplitude T(r, x) with decreasing x 
is described by the B-JIMWLK equation ll2 , 
established in the leading ln(l/x) approxi- 
mation, and by its mean-field version the BK 
equation 3 . 

What the dipole sees is what is pictured 
in Fig.l with kx — > 1/r. In the leading- 
twist regime T <C 1, and in the saturation 
regime T = 1. From the BK equation and 
its equivalence with the FKPP equation, well 
known in statistical physics, it was shown 
that the growth of the dipole amplitude to- 
wards the saturation regime occurs in a par- 
ticular way 4 . When approaching the satu- 
ration regime, instead of being a function of 
a priori the two variables r and x, the dipole 
scattering amplitude is actually a function 
of the single variable rQ s (x) up to inverse 
dipole sizes significantly larger than the sat- 
uration scale Q s {x). 

This means that in the diagram shown 
Fig. 1 , lines parallel to the saturation line are 
lines of constant densities, along which T is 
constant. This property which extends to- 
wards small values of r (or large values of Q 2 ) 
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Fig. 2. The cross section crJ f^ X as a function of 
T = Q 2 /Q 2 a {x) for x<0.01. 

manifest itself in the data. Indeed, this im- 
plies the geometric scaling of the total cross- 
section at small x : 

air x (x,Q 2 ) = *ir x (T=Q 2 /Q 2 s (x)) . 

(i) 

As shown in Fig. 2, this has been con- 
firmed by experimental data 5 with Q s (x) — 
(x/xo)~ x / 2 GeV and the parameters A = 
0.288 and x = 3.04 1CT 4 taken from 6 . 

The scaling law T(r, x) = T(rQ s (x)) also 
implies the geometric scaling of the diffrac- 
tive cross section at small xp and fixed (3, 
with the same saturation scale: 

ajlff Xp ((3,r d = Q 2 /Q 2 s (x P )).(2) 

This is in agreement 7 with the data from 
HERA (see Fig. 3). More quantitative anal- 
ysis have been carried out: saturation mod- 
els 6 ' 8 ' 9 fit well cr^ot^ data and give accu- 
rate predictions for aJ i j'p Xp , and more ex- 
clusive processes like deeply virtual Compton 
scattering and vector meson production 10 . 

More studies contributed to the success 
of saturation in the context of forward par- 
ticle production in hadron-hadron collisions. 
When producing particles (jets, pions, dilep- 
tons, heavy flavors) with transverse momen- 
tum pt and rapidity y, the typical values of 
x probed in the hadrons wavefunctions are 
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Fig. 3. The cross section do^Tf* p /d/3 as a func- 
tion of Td = Q 2 /Qs(a;p) in bins of /3 for xp<0.01. 

xi,2~Pt e ±y / v /s. Therefore particle produc- 
tion at forward rapidities and high energies 
involves high values of x for one hadron and 
small values of x for the other. To describe 
such a process, the dipole picture can also be 
used 11 , allowing many investigations 12 . 

Among other things, saturation pre- 
dicted the high-rapidity suppression of the 
nuclear modification factor in d-Au colli- 
sions, which was observed at RHIC. Other 
predictions, such as the suppression of back- 
to-back correlations, are in qualitative agree- 
ment with the data. 

Recently, some limitations of the B- 
JIMWLK equation were pointed out 13 ' 14 . 
This triggered a series of papers that aimed 
at completing the B-JIMWLK equation by 
including Pomeron loops 15,16 in the high- 
energy evolution. These studies improved 
our understanding of different aspects of 
high-energy scattering. For instance, they 
allowed the constructions of an effective ac- 
tion 17 , a generalized dipole model 18 and a 
reggeon field theory 19 . 

In the following, we discuss the phe- 
nomcnological consequences of Pomeron 



loops. Within the large— N c limit, the dipole 
amplitude T(r, x) can be obtained from a 
Langevin equation 15 that features a noise 
term, similar to the stochastic FKPP equa- 
tion 14 . Its solution is an event-by-event 
dipole scattering amplitude T(rQ a (x)) char- 
acterized by a saturation scale Q s which is 
a stochastic variable. Using results derived 
in the weak-noise 20 and strong-noise 21 lim- 
its, it was shown that ln(Q s /Qo) (with Qq a 
reference scale) is distributed according to a 
Gaussian probability law 22 . Denoting Y = 
ln(l/a;), the average value is \n{Q 2 s / Q\\) = \Y 
and the variance is tr 2 = DY. 

One obtains the physical amplitude T 
after averaging T over the different real- 
izations of the noise. The dispersion co- 
efficient D allows to distinguish between 
two energy regimes: the geometric scaling 
regime (DY<^1) and diffusive scaling regime 
(D7> 1). In formulae, the averaged ampli- 
tude features the following scaling behaviors: 



T(r,Y) 



Y<Sil/D 



,00) . 



(3) 



T(r,Y) Y>> = /D T ds (]n(rQ s (Y))/VDY) .(4) 

In the saturation region rQ s > 1, T(r, Y) = 1. 
As the dipole size r decreases, T(r,Y) de- 
creases towards the dilute regime, following 
the scaling laws (3) or (4) , depending on the 
value of DY. 

In other words, Pomeron loops turn the 
geometric scaling regime into an intermedi- 
ate energy regime, in which the dispersion 
of the events is negligible. It predicts the 
emergence of a new scaling law at very high 
energies, when the event dispersion is impor- 
tant: diffusive scaling 23 . This is summa- 
rized in Fig. 4 (with r — + 1/Q, Fig. 4 applies 
also to DIS). A key feature of the diffusive 
scaling regime is that, up to inverse dipole 
sizes much larger than the average saturation 
scale, cross-sections are dominated by events 
which feature the hardest fluctuation of the 
saturation scale. In average the scattering is 
weak yet saturation is the relevant physics. 
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Fig. 4. A phase diagram for the high-energy limit of 
DIS in QCD. Shown are the average saturation line 
and the boundaries of the scaling regions at large val- 
ues of ln(l/r 2 ). With increasing Y, there is a gradual 
transition from geometric scaling at intermediate en- 
ergies to diffusive scaling at very high energies. 

While it seems that HERA and RHIC 
are probing the geometric scaling regime, 
this has important implications for the LHC. 
In the context of forward particle produc- 
tion 24 , the energy there might be high 
enough to reach the diffusive scaling regime. 
However, our poor theoretical knowledge of 
the dispersion coefficient D prevents more 
quantitative statements at the moment. 
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